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Abstract 

Expressions for dependences of the pre-exponential factor 0-3 and the thermal activation energy £3 of hopping electric conductivity of holes 
via boron atoms on the boron atom concentration N and the compensation ratio K are obtained in the quasiclassical approximation. It is 
assumed that the acceptors (boron atoms) in charge states (0) and (—1) and the donors that compensate them in the charge state (+1) form a 
nonstoichiometric simple cubic lattice with translational period Rh = [(1 + K)N]~ 1 '' A within the crystalline matrix. A hopping event occurs 
only over the distance Rh at a thermally activated accidental coincidence of the acceptor levels in charge states (0) and (—1). Donors block 
the fraction K/(l — K) of impurity lattice sites. The hole hopping conductivity is averaged over all possible orientations of the lattice with 
respect to the external electric field direction. It is supposed that an acceptor band is formed by Gaussian fluctuations of the potential energy 
of boron atoms in charge state (—1) due to Coulomb interaction only between the ions at distance R\,. The shift of the acceptor band towards 
the top of the valence band with increasing N due to screening (in the Debye-Hiickel approximation) of the impurity ions by holes hopping 
via acceptor states was taken into account. The calculated values of a-i(N) and E3(N) for K ~ 0.25 agree well with known experimental data 
at the insulator side of the insulator-metal phase transition. The calculation is carried out at a temperature two times lower than the transition 
temperature from hole transport in u-band of diamond to hopping conductance via boron atoms. 

Key words: A. Boron doped diamond, B. Moderate compensation, C. DC hopping conductivity, D. Activation energy, E. Nearest neighbor hopping 
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1. Introduction 

Following the observation of helium temperature supercon- 
ductivity in heavily boron doped diamond, numerous studies 
of this material were performed (see, e.g., reviews [1-3]). The 
possible use of intermediately boron doped diamond in semi- 
conductor applications [4] justifies studies of its conductivity at 
room temperature. With temperature lowering from the room 
temperature to the temperature of liquid nitrogen a conduc- 
tion of holes in u-band (propagating regime [5]) changes into 
hopping conduction of holes via boron atoms. The dc hopping 
conduction in boron-doped diamond is observed in the dark at 
significantly higher temperatures and concentrations of boron 
than in silicon and germanium crystals doped with the same 
acceptor impurity at comparable compensation ratios. Because 
of the progress in the synthesis technology of high-quality ho- 
moepitaxial crystalline diamond films with controllable boron 
doping, reliable experimental data on the hopping conductiv- 
ity Oh via boron atoms comparable with computational models 
was eventually obtained [6-8]. However, a satisfactory quanti- 
tative description of the hopping transport of holes in diamond 
crystals is still lacking. 

In this paper we limit our consideration to the hole hopping 
regime via nearest neighbor boron atoms (NNH regime). In this 
case, the dc hopping conductivity is (see, e.g., Ref. [9]) 



er h = cr 3 exp 



£3 



(1) 



where 03 = I//03 is the pre-exponential factor, e^/k^T is the 
ratio of the activation energy of hole transport via impurity 
atoms to the thermal energy. 

When describing small polaron hopping over the lattice sites 
in an ionic crystal, Holstein (see reprint [10]) introduced the 
concept of a "coincidence event" for polaron potential wells. 
In the model [10] (see also Refs. [11,12]), a polaron hop is 
assumed to occur when the energies of the initial occupied state 
and the final vacant state coincide. 

In Ref. [13], a model of fluctuation-induced "alignment" of 
the energy levels of localized states (of impurity atoms) due to 
electron-electron interaction was proposed to describe the dc 
conduction of doped semiconductors. It was assumed that tem- 
poral fluctuations of the energy of localized states are caused by 
hopping diffusion of electrons via these states. Another model 
(the variable range hopping (VRH) conduction model [14]) has 
been proposed, assuming that electron (or hole) hops occur via 
resonance tunneling between atoms of majority impurity. En- 
ergy levels of two impurity atoms enter into resonance due to 
Coulomb potential fluctuations induced by stochastic changes 
in the occupation state of other doping impurity atoms. l How- 
ever, in studies [13,14] numerical calculations of the experi- 
mentally observed quantities were not carried out. 
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1 The model of the fluctuation-induced preparation of a barrier through which 
an atom (or even a molecule) can tunnel made it possible to explain the main 
characteristics of solid-phase cryochemical reactions (see, e.g., [15]). 
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In Refs. [16,17], the dc hopping conductivity and thermo- 
electric power in the NNH regime were described for neutron 
transmutation doped p type germanium. It was supposed that 
the doping impurity (acceptors) and a compensating impurity 
(donors) form a common "impurity lattice" in the crystalline 
matrix. It was assumed that hole hopping occurs only during 
thermal-induced alignment (coincidence) of energy levels of 
the acceptors in charge states (0) and (—1). Donors, which are 
all in charge state (+1), block some sites of the impurity lattice. 
At the instant of energy level coincidence (up to a broadening 
of the acceptor energy levels due to the finite time of hole local- 
ization on the acceptor) of a neutral and a negatively charged 
acceptor, a "resonant" two-site cluster is formed: the hole on 
acceptor 1 becomes bound to a negatively charged acceptor 2 
and belongs simultaneously to these two acceptors. After some 
time, the resonance conditions are no longer satisfied and the 
hole can become localized on acceptor 2 or remain on acceptor 

1 . After that, the acceptors 1 and 2 can again form a resonant 
cluster or form resonant clusters (acceptor pairs) with other ac- 
ceptors. 

The aim of our work is to determine the dependences of 03 
and £3 in Eq. (1) on the concentration of boron atoms in mod- 
erately compensated, intermediately doped diamond crystals. 
For this purpose the model from Ref. [16] is developed. 

2. Statistics of hydrogen-like impurities 

Let us consider a p type uniform crystalline semiconductor 
with acceptor concentration N = No + N_\, where No and 
A_i are the concentrations of acceptors in charge states (0) 
and (—1), respectively. There are also donors, all in charge 
state (+1), with concentration KN, where < K < 1 is 
a compensation ratio of acceptors by donors. 2 The electrical 
neutrality condition has the form 



KN = N_!. 



(2) 



According to Eq. (2), the probability that a randomly chosen 
acceptor is in charge state (—1) or (0) is equal to K or (1 — K), 
respectively. Thus, the average (over the crystal) fraction of 
acceptors in charge state (—1) is [18-20] 
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where Q is the distribution of acceptor energy levels E with 
respect to the average value E when /_i is the probability that 
acceptor with energy level E is ionized. Further, we assume 
that Q is a Gaussian distribution: 
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where W is the effective width of the acceptor band. According 
to Refs. [18,21] the probability that an acceptor with energy 



2 We consider the boron concentrations N much less than the concentration 
at which in diamond the transition from insulator state to metallic one (Mott 
transition) TVm rj 2-10 20 cm' 3 is observed [6]. The calculations of Nu 
dependence on compensation ratio K were carried out in the work [18]. 



level E > above the top of the w-band (E v = 0) of an 
undoped crystal is ionized can be written as 
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here f3 = 6 is the degeneracy factor of the energy level of a 
boron atom in diamond, Ep is the Fermi level relative to the 
top of the u-band (Ep < in the band gap of diamond), k%T 
is the thermal energy, u = (E — E)/k B T and £ — ln/3 = 
(Ep + E)/ksT are dimensionless acceptor energy level and 
Fermi level relative to the center E of the acceptor band. 

We assume that the doping impurity (hydrogen-like accep- 
tors) and the compensating impurity (hydrogen-like donors) 
form a common nonstoichiometric simple cubic lattice within 
the crystalline matrix (cf. [22]). The translational period of this 
lattice is R h = [(1 + i^Q/V] -1 / 3 , where (1 + K)N is the total 
concentration of impurities. In further consideration we sup- 
pose that Rh is the length of hole hop between acceptors in 
charge states (0) and (—1) in the impurity lattice. 

The acceptor band width, taking into account the Coulomb 
interaction of the ionized acceptor with the ions in the first 
coordination sphere of the impurity lattice, is equal to [16,17] 

"-(rw?) = i^ ITT? • <« 

X 2 — 1 / X ' 



where Vi = 2K/(1 + K) = 2KE is the probability that any 
of 6 sites of a simple cubic impurity lattice (in the first co- 
ordination sphere) near the z-th ion is occupied by an ionized 
acceptor or donor, 5 = 1/(1 + K) is the correlation factor, 
i.e., the fraction of majority impurity at impurity lattice sites, 



Pi\ 



/(AiteRh) is the magnitude of the Coulomb interac- 



tion energy between the i-th ion and an ion in the impurity 
lattice at distance Rh from it, e is the modulus of the electron 
charge, s = 5.7£o is the static permittivity of diamond, £0 is 
the electric constant. For derivation of Eq. (6) it is taken into 
account that the Coulomb interaction energy of the ions in the 
first coordination sphere averaged over the impurity lattice is 
equal to zero: ^2 i=1 Villi = 0. 

According to Eq. (6), we assume that the acceptor band is a 
"classical" one, i.e., the spread of the energy levels of the boron 
atoms is much greater than the quantum-mechanical broadening 
of these levels due to the finite time of hole localization on the 
acceptor. 

The position of the acceptor band center E relative to the 
top of the w-band (E v = 0), according to Refs. [18,19], is 

__ 3e 2 

E ~ I ~ 167T£(A h + R h ) ' (7) 

where / = 370 meV is the energy level of a single accep- 
tor (boron atom in diamond), Ah is a screening radius of the 
Coulomb field of the ion, R h = [(1 + A^A]- 1 / 3 is the min- 
imal possible distance between ions in the impurity lattice. In 
the Debye-Hiickel approximation, according to Refs. [18-20], 
we find 



Ar 



l K{l-K)N 



(8) 



£k B T£,h ' 

where ^h ^ 1. In accordance with Ref. [20], the reciprocal value 
of the quantity §, in the modified (generalized) Einstein relation 



for hopping migration of holes via acceptors is determined by 
expression: 
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where Dh/Mh is the ratio of the diffusion coefficient of hopping 
holes to their mobility. 

For the narrow acceptor band (W <C k B T), taking into ac- 
count Eqs. (4)-(6) it follows from Eq. (9) that the Einstein rela- 
tion is fulfilled in the classical form (i.e., £h = 1)- In this case, 
we find from Eq. (2) the quantity < = In j3 + (E F + E)/k B T w 
— ln[K/(l — K)] taking into account Eqs. (3)-(5). 

For the wide acceptor band (W 3> k B T), according to 
Ref. [20], we obtain from Eq. (9): 



Let us take into account all possible orientations of the cubic 
impurity lattice with respect to the direction of the external elec- 
tric field S '. The space of possible orientations is a semi-sphere 
with a normalized element of the surface (1/2tt) s'm9d(pd9. 
All edges of the impurity lattice unit cells, each of length R^, 
directed at angle 6 to the field strength § make a contribution 
<7hcc cos 9 to the conductivity, where a^ x is given by Eq. (12). 
As a result, taking into account Eq. (1), we find: 

r2iv fw/2 
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cos 9 sin 9d9 = 
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(13) 



£h « K{1 - i.)7V2^cxp(C727 2 ), 

where 7 = W/k B T ^> 1, In this case, Eq. (2) takes the form: 
2K « 1 - erf(C/\/27). 



3. Hopping current and conductivity 

Let us consider the range of N, K and T values at which the 
NNH regime is only realized, and £3 weakly decreases upon a 
reduction in temperature. 

According to Refs. [16,17,23,24], a stationary hopping cur- 
rent density Jh of holes over acceptors in a sample subjected to 
an external electric field of the strength $ = —dip/dx directed 
along the x axis has the form: 

dN 



./„ = < _Y h ( Mh^-TV^- In -^2- 
dx iv_i 



To obtain the dependence of hopping conductivity <7h on the 
values of JV, K and T using Eqs. (13) and (12) one needs to 
find the frequency of hole hopping Th in the impurity lattice. 

For the not very small and not very large ratios of acceptor 
compensation by donors (tentatively for 0.1 < K < 0.9) the 
correlation between the position of an acceptor in the impurity 
lattice with translational period R^ and the acceptor energy 
level E may be neglected. Following Ref. [16], we suppose 
that every hole hop of length R^ between acceptors 1 and 2 
in charge states (0) and (—1) occurs only when their energy 
levels (Ei =E + uik B T and E 2 = E + u 2 k B T) accidentally 
coincide. The number of hole transitions between acceptors 
(boron atoms) per energy level coincidence event is equal to 
the integer part of the ratio of duration tu (u) of an event of 
coincidence of levels (m = u 2 = u = (E T — E)/k B T) to the 
time t(u) of a tunneling event. We assume also that for the time 
interval t the total duration of all events of level coincidence 



where N^ = NqN-i/N is the effective concentration of holes 
which hop between acceptors in charge states (0) and (—1), 
Mh = S_R^ dT/dip > is the hopping mobility of holes, 
S = 1/(1 + K) is the fraction of acceptors at impurity lat- 
tice sites, R b = [(1 + K)N]~ 1/3 is the length of the hole hop, 
(— dr/dx)i?h is the difference of the average hole hopping fre- 
quency in the direction along and against the external electric 
field, (p(x) is the electric potential, D^ = Si?^Th/6 is the dif- 
fusion coefficient, Th/6 is the average hole hopping frequency 
along one of the six directions (along the edges) of the impu- 
rity lattice for zero external field. The hopping conductivity ahx 
of holes in the impurity lattice for external field S orientation 
along the edge of the unit cell of the cubic impurity lattice is 
given by Eq. (10), taking into account Eq. (2), in the form 



= Uhx^—eDh— — , (10) i s equal to t c (u) = ^2 k tk(u). We approximate the conditional 

probability that exactly j transitions of a hole occur at the 
coincidence of levels of two nearest neighbor acceptors by the 
Poisson distribution [25,26]: 

t c (u)~ 
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where t c (u)/r(u) = ^2^L jP{j\u} is the average number of 
hole transitions between nearest acceptors (boron atoms), t(u) 
is the duration of a tunneling transition of a hole from the 
neutral acceptor to the ionized one, j = 0, 1, 2, 3, . . .. 

Thus, the frequency of hole hopping between two acceptors 
at accidental alignment of their energy levels E T = E + uk B T 
(an average number of hole transitions for the time t) is 
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Let us find from Eq. (9), using D^ = Si?^Th/6, the hole mo- 
bility Mh in terms of the equilibrium frequency Th of their 
hopping via acceptors 

Mh = -^l. = eSJ? h 2 r h 

S,hk B T 6£hk B T ' 

Thus, from Eq. (11) we obtain the hopping conductivity due 
to hole hopping with hop length R^ along the external electric 
field (along the x axis) 

e 2 EK(l - K)NRlT h 



From the theory of Markov chains [25,26], it follows that, if 
the hole transitions between two acceptors are observed over a 
long time interval (t ^> r(w)), then the fraction of time spent by 
the acceptors in one of two possible states (when their energy 
levels are coincident or noncoincident) is approximately equal 
to the stationary probability of the acceptors being in these 
states. Thus, the ratio t c (u)/t <C 1 is approximately equal to 
the probability that the energy levels of two nearest neighbor 
acceptors in charge states (0) and (—1) are aligned (energy 
values belong to interval (u, u + du)): 
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where the product f (w)/-i (u) is obtained using Eq. (5) and the 
distribution density G(u) is given by Eq. (4) at the substitution 
of E — E by uk B T. The quantity P(u) du in formula (16) gives 
the conditional probability that energy levels of a randomly 
chosen pair of acceptors in the charge states (0) and (—1) belong 
to the interval (u, u + du) (see Appendix). 

Further, we take into account that the energy level E T (u) — 
E + uk B T > above the top of the w-band is associated 
with a radius a T (u) = e 2 /(8ttsE t ) of hole localization on the 
acceptor with ionization energy E T . The Bohr radius for the 
center (u = 0) of the acceptor band is o T (0) = e 2 /(8iTeE). 

Within the framework of the theory of the hydrogen molec- 
ular ion (Hj) [27] a duration of hole tunneling between two 
acceptors at the distance i?h when their energy levels coincide 
(u\ — U2 — u) can be estimated as [28,29]: 

T ™ = JE^u-y (17) 

where h is the reduced Planck constant, 8E T (it) is the broad- 
ening (splitting) of the energy levels E T (u) = E + uk B T = 
e 2 /[87rea T (w)] of the acceptors when the hole tunnels between 
them: 



SE T (u) = AE T (u) 
^ p(l + p) cxp(-p) - [1 - (1 + p) cxp(-2p)]S 



(18) 



p(l S 2 ) 

p(u) = R h /a T (u), S(u) = [l + p+ (p 2 /3)] exp(-p). 

Let us average T(u) over the distribution of energy levels 
which form the acceptor band with effective width W (see 
Eq. (6)). Taking into account Eqs. (15) and (16), the average 
hole hopping frequency Th between two acceptors at the dis- 
tance i?h in the impurity lattice can be written in the form: 

+oc f +oo p( u ) 
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(19) 



The duration of hole tunneling t(u) according to Eqs. (17) 
and (18) monotonically increases when the tunneling level 
E T moves deeper into the band gap. The function P(u) ex 
G(u)fo(u)f-i(u) has its sharp maximum at value u = u m 
which satisfies the equation 



+ 7 2 tanh[(u m + 0/2] = 0, 



(20) 



where 7 = W/k B T, C = In j3 + (E F + E)/k B T. This allows us 
to take t(u) out of the integral in Eq. (19) at u = u m denoting 
r(u m ) = 73. It follows from Eq. (20) that if temperature T —> 
then Um — *■ — C> an d if acceptor band width W — > then 
w m — *• 0. Thus, taking into account Eq. (9), the average (over 
the impurity lattice) frequency of hole hopping is 
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(21) 



where T^ — l/r(u m ) = I/73 is the frequency of hole tunnel- 
ing between acceptors in charge states (0) and (—1) with en- 
ergy levels E T (u m ) = E + u m k B T. The electrical neutrality 
condition (2) can be solved for the dimensionless Fermi level 
( - In j3 = (E F + E)/k B T. Calculating the position of the ac- 
ceptor band center E using Eq. (7), we find the Fermi level 



value E F = k B T(( - ln/3) - E < relative to the top of the 
w-band. The value u m is calculated from Eq. (20). 

We find the pre-exponential factor in the temperature depen- 
dence (1) from expression (13) taking into account Eqs. (12) 
and (21): 
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(22) 



where the duration of hole tunneling I/T3 = T3 = r(w m ) is 
determined using Eq. (17) for u = u m from Eq. (20), and the 
factor £h ^ 1 is given by formula (9). 

The activation energy of hopping conductivity follows from 
Eq. (21) using Eqs. (12) and (13) in the form: 



£ 3 = -fc B Tln(r h /r 3 ) = k B T\n& 



(23) 



It follows from Eq. (23) using Eq. (9) that £3 — *■ at the 
acceptor band width W — > . 

According to Eq. (23) and taking into account Eq. (9), the 
activation energy £3 decreases sublinearly upon lowering the 
temperature T, because £h increases in this case. Let us define 
the characteristic temperature value 2], at which £3 is measured. 
For some temperature Ij the conductivity a p of holes in the v- 
band is equal to the hopping conductivity Ch of holes via boron 
atoms. The dependence of the temperature Ij on the concen- 
tration N of boron atoms is determined in the same way as in 
Ref. [18], i.e., from the dependences of the logarithm of the 
total electric conductivity a = <r p + Ch on the reciprocal tem- 
perature 1/T for different A^ and 0.05 < K < 0.5, using ex- 
perimental data [6-8]. 3 We derived the numerical dependence 



of T\ on N in the form 
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(24) 



where [T-] = K, [N] = cm" 3 . 

The average temperature 2],, at which the values of 0-3 and 
£3 were measured in the experiments [6-8], can be estimated 
as Th = Tj/2, i.e., it is assumed to be equal to the average 
value in the interval from absolute zero to Tj. (It is likely that in 
the temperature interval from Tj to 7], the NNH regime of hole 
hopping via boron atoms dominates, and for the temperature 
below Th the VRH regime occurs.) 

In Fig. 1, the dependence of the reciprocal value of the pre- 
exponential factor (22) for hopping resistivity (^3 = 1/173) 
on the concentration N of boron atoms in diamond (for the 
compensation ratio K = 0.25 and temperature Th = Tj/2) 
is shown by the solid line. When we change A -1 / 3 from 2 
to 11 nm the temperature Th changes, according to Eq. (24), 
between 160 and 90 K. In this case the ratio W jk B Th changes 
approximately from 15 to 5. Thus, according to Eq. (20) we 
obtain u m s=s — (. 

It should be noted that the calculation according to Eq. (22) 
also agrees well (see Ref. [16]) with the data for p$ in/? Ge:Ga 
at K = 0.35 and T h = Tj/2, where according to Ref. [19], the 



3 In the works [6,7] the total concentration N = A^o + N—\ of boron 
atoms in diamond was determined from the secondary ion mass spectroscopy 
measurements. In the work [8] the concentration Nq of electrically neutral 
boron atoms was measured from the one-phonon band of IR absorption at 
T = 300 K. Then we estimated for K f» 0.25 the boron concentration to 
be N = N /(l - K). 
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Fig. 1. Dependence of the pre-exponential factor P3 = 1/0-3 on me concen- 
tration N of boron atoms in diamond. Solid line is our calculation using Eq. 
(22) for the compensation ratio K = 0.25 at the temperature Tj, = Tj/2 
defined by Eq. (24). Experimental data of papers [6], [7], [8] are denoted by 
a, b, c, respectively 

temperature Tj = 5.3-10 4 7V 027 , [Tj] = K, [N] = cm" 3 . How- 
ever, in this case, the acceptor energy level of the single gal- 
lium atom in germanium is I = 11.3 meV, relative permittivity 
is e/eq = 15.4, and the value N^ 1 ^ 3 changes approximately 
from 30 to 140 nm. 

In Fig. 2, the dependence of the activation energy £3 of 
hopping conductivity on the concentration N of boron atoms 
in diamond for compensation ratio K = 0.25 is shown by 
the solid line. The calculation was performed using Eq. (23) 
taking into account Eq. (9) at the temperature TJ, = Tj/2. In 
Fig. 2, the calculation of £3 « 0.7e 2 A fl / 3 /(47T£) according 
to the model [9], for K ss 0.25 and relative permittivity of 
diamond e/eq = 5.7, is shown by the dashed line. It can be 
seen that the calculation according to Eq. (23) agrees well with 
the experimental data, while the calculation using the model [9] 
gives overestimated values of the activation energy £3. 

Note that, for all covalent semiconductors with hydrogen-like 
impurities, with the increase of their concentration, the quantity 
£3 in Eq. (1) decreases after reaching a maximum value (see, 
e.g., Ref. [9]). The model [16] (see also outlines [30]) allows 
one to describe such a decrease at the expense of the broadening 
8 E T of the energy levels of single acceptors with the increase 
of their concentration due to the finite time of localization of 
hole on the acceptor. In this case, the realization of the acceptor 
level resonance condition (16) becomes more likely, because 
8E T defined by Eq. (18) becomes comparable with acceptor 
band width W given by Eq. (6). 

4. Conclusions 

We have developed a model of dc hopping conductivity of 
holes via acceptors (boron atoms) in diamond crystals. It was 
supposed that boron atoms with concentration N and donors 
with concentration KN form a common nonstoichiometric 
simple cubic lattice within the crystalline matrix. The transla- 
tional period of impurity lattice is Rh = [(1 + -ft^iV] -1 / 3 . In 
this case, only hole hops of length Rh between acceptors at the 
instances of accidental alignment of their energy levels due to 
the thermal fluctuations were taken into account. The hopping 
conductivity was averaged over all possible orientations of the 
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Fig. 2. Solid line is the dependence of the activation energy £3 of hopping 
conductivity on the concentration N of boron atoms in diamond, calculated 
by Eq. (23) for the compensation ratio K = 0.25 at temperature Tj, = Tj/2 
defined by Eq. (24). Experimental data of papers [6], [7], [8] are denoted by 
a, b, c, respectively. Dashed line is the calculation using the model [9] 

cubic impurity lattice with respect to the direction of the exter- 
nal electric field. It was taken into account that the energy levels 
of boron atoms are distributed due to the Coulomb interaction 
between the ions in the first coordination sphere of the impu- 
rity lattice. We assumed that the acceptor band width is much 
larger than the quantum broadening of acceptor energy levels 
because of the finite time of hole localization on boron atoms. 
Our calculations of the pre-exponential factor 0-3 and the acti- 
vation energy £3 of the hole hopping transport o^ over boron 
atoms depending on their concentration agree well with exper- 
imental data [6-8] for moderately compensated (K « 0.25) di- 
amond crystals. The calculation was performed for temperature 
Ti = Tj/2, two times lower than the transition temperature Tj 
from w-band dc hole conduction regime to hopping over boron 
atoms. 
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Appendix 

To obtain formula ( 1 6), let us supplement the notations of the 
main text of the paper. The probability of finding an acceptor in 
charge state (—1) is equal to P{ — 1} = K, and the probability 
of finding an acceptor in charge state (0) is equal to P{0} = 
1 — K. The conditional probability that an acceptor with energy 
level s = u is in charge state (0) is equal to T{0 \s = u} = 
fo(u). The conditional probability that an acceptor with energy 
level s = u is in charge state (—1) is equal to P{ — 1 | S = u) = 

The conditional probability of an event A occurring, provided 
that event B has occurred, has a form [25,26]: P{yl|_B} = 
P{AC\B}/P{B}. Hence, exchanging A and B, we have P{An 
B} = P{B\A}P{A}, and consequently 



P{A\B} 



P{B\A}P{A} 
PjB} " 



(A.l) 



Formula (16) has a form (A.l), where the event B corre- 
sponds to one acceptor in an acceptor pair being in charge state 
(0) and the other being in charge state (—1); the event A corre- 
sponds to each acceptor in the pair having the energy level u. 
Thus, 

P{B} = P{pair(0, -1)} = 2K(1 - K) (A.2) 

and the conditional probability 

P{B\A} = P{pair(0, -1) | s = u} 

= 2P{{) | a = u}P{-l \s = u} = 2/ (u)/_i(«). (A3) 

Because the distribution of acceptor energy levels relative 
to their average value (s = 0) is continuous, instead of the 
probability -P{A} we use the probability of finding the acceptor 
energy level s within the interval (u, u + du): 

P{se{u,u + du)} = g(u)du. (A.4) 

Now the conditional probability to find a pair of acceptors in 
charge states (0) and (—1) with energy levels within the interval 
(u, u + du) can be obtained by substitution of Eqs. (A.2)-(A.4) 
into Eq. (A.l): 

20(u)/o(u)/-i(u)du 



P{s e {u,u + du) | pair(0, -1)} = 



2K(1-K) 



Thus, the quantity 



/o («)/-i («) 



Q(u)du = P(u)du gives 



K(l-K) 

the conditional probability that a pair of acceptors in charge 
states (0) and (—1) have their energy levels within the interval 
(u, u + du), i.e., we obtain formula (16). 
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